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I . INTRODUCTION 


Weather radar installations have generally been confined to 
frequencies at which the attenuation of the signal by liquid water and 
the atmosphere is not very significant. The difficulties of making proper 
attenuation corrections and the instability in current methods, which can 
result from errors in calibration, have discouraged their use at attenu- 
• ating frequencies D,2]. Precise calibration of radars has long been 
recognized as a major problem. Traditional calibration n^thods have 
been based upon the use of a known target, such as a sphere [3], or 
comparison with ground rain rate measurements [4]. 

The experiment described herein uses three radars, two of which 
operate at attenuating frequencies. The three-frequency High Resolution 
Radar/Radiometer System at the Ohio State University includes radars in 
the S, X and Ku bands, (3.05, 9.35 and 15.56 GHz), looking simultaneously 
at the same spatial region. Such a configuration possesses son« unique 
advantages. 

This system was developed to explore the potential improvement in 
radar rain rate measurement accuracy which would result from calibration 
of radar backscatter at a second frequency or in terms of an indepen .-nt 
measurement of radiometric temperature. The use of simultaneous measure- 
ments relaxes the need for high absolute accuracy in each individual 
instrimient as well as making possible compensation for variations in 
drop-size distribution. The use of higher frequencies also improves the 
spatial resolution capability of the radar. 

The Ohio State University Multi -Frequency High Resolution Radar/ 
Radiometer System was implemented and operated for a brief period of 
time during the calibration and debugging phase of the program in 
1973. Subsequently, a shift in program emphasis and funding curtailed 
this unique effort. Nevertheless, during this "shake-down" period, 
approximately 100,000 records of high resolution radar/ radiometer data 
were recorded. Editing and analysis of these data have been performed 
and are described herein. 

In order to accomplish this objective, several calibration tech- 
niques were examined and tested. The calibration procedure finally 
chosen is also described in the following. Attenuation corrections were 
applied and the resulting multi -frequency data were checked for self 
consistency. The techniques and theory used follow some new theoretical 
work on weather radar calibration [6]. 

Cell width, equivalent reflectivity and peak equivalent reflectivity 
distributions are computed using the edited and calibrated data. The 
average cell width at several relative reflectivity levels below the 
peak are obtained for the cells observed. The feasibility of using 
attenuating frequency radars is demonstrated. 
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II. EQUIPMENT AND FACILITIES 


Tte High Resolution Radar/Radi cpeter System at The Ohio State 
University consists of three radars operating at 3.05 GHz (S-Band), 

9.35 GHz (X-Band) and 15.56 GHz (Ku-Band), together with radiometers at 
3.30, 8.50 and 15.10 G^z respectively. The radars are fired in sequence 
under computer control, the pulse-repetition frequency of each radar • 
being 100 Hz. The initial range is adjustable from 15 to 45 km and the 
range window is 30 km with a resolution of 300 meters. The returns of 
32 pulses of each radar are sunned on a range bin-by-range bin basis 
and recorded on digital magnetic tape, together with time, look angle 
and other relevant information. A summary of the system parameters is 
given in Table 1. For a complete Ctescription see reference [5]. 


III. EXPERIMENTAL CONDITIONS 

The available data consist of about 100,0(K) records. These were 
taken during the final stages of construction and commissioning of 
the equipment while debugging and calibration were being performed. 
Consequently, hardware problems are evident in some of the data. Further- 
more, the termination of the effort curtailed the calibration and 
documentation process. Some uncertainty remains in parameters such as 
antenna gain, transmitted power and receiver gain slopes. 

The problem of data analysis is further aggravated by the four 
year lag between the taking of the data and its analysis. Some 
personnel working on the system, who might have been able to provide 
valuable Information on the problems encountered, are no longer 
available. Nevertheless, it was felt that the information contained in 
the available data warranted its analysis and that adequate calibration 
information was available to make a complete calibration possible. 


IV. THE DATA 

Approximately 100,000 records on magnetic tape were edited to 
remove unsuitable data. Records showing hardware problems such as loss 
of transmitter power, integrator overflow, etc. were discarded. Many 
records showing only portions of cells of significant size were also 
rejected. The remaining usable 10,000 records show four rainfall events. 
The X- and Ku-faand radars were operational during all four data periods, 
while the S-band radar was operated during only one of these. All three 
radars were oriented to the same pointing angle, which was fixed in 
elevation at 6®,during each event. The radars were allowed to dwell 
on a rain cell for some time and then rotated in azimuth to another 
cell of interest. Seven rain cells are seen as they move through the 
fixed beams. 
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V. 


SYSTEM REPRESENTATION 


A block diagram of one of the receivers Is shown In Figure 1, 

Based on available preliminary calibrations each radar receiver was 
represented by a straight line characteristic on a log-power In vs. 
output voltage plot with gradient GR dB/volt and noise level FN^ef dBm. 
After the A/D converter and 32 pulse digital summer, the overall gradient 
was A dB/dIgItal level recorded on tape. If 

= received power per pulse from any range bln, 

D = digital level for that range bln, 

D = digital level at the tall edge of record, representing 
° receiver noise, 

then It can be shown that 


(iTW) 


= FN 


ref 


(mw) 


(D-Do^fo 
10 ° - 1 


0 ) 


Also, the monostatic beam-filled radar equation can be written as 




-2A_ 


■eq 


where 


P. CT 

Q = 

1024 An2 




r? + 2 
c 


P^ = transmitted power 
G = antenna gain 
0 = antenna beam width 

T = pulse width 
A = wavelength 
c = velocity of light 

n„ = complex refractive index of water 

^ 6 3 

Zeq = equivalent reflectivity factor mm /m 

r = range to point of reflection 

A^ = total one way attenuation to point of reflection. 


( 2 ) 


(3) 
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Block diagram of a radar receiver. 






It has been shown by Hodge [6] that If Zgq and the specific attenuation, 
a (dB/toi), could be represented In terms of rain rate, R (tm/hr), as 

b. 


O = R 


'RD 


(5) 


then 


-2A 


e ^ = (1 - I^) 


"z/^RD 


where 


I « 
r 


^ ®RD f / 2p v^RD''*^Z . 

4.343(a,Q) ^ "o 


Hence 




“z, 


eq 


( 6 ) 


(7) 


( 8 ) 


Substituting from Equation (1) for P^, dividing both sides by FN^ef 
and writing 

2 


“ "7 

f\ef 


2 2 3 

P*G^e^CT n'’ 


X'^FNref 


«2 1 
"c * 2 


(9) 
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, 9^ - I,)'''''” Z,, 


2 a F 

. ^ *RD ^ f 

P b„/b" 1 

4.343(a2Q’) ^ o 


( 10 ) 




-| I^RD^’^Z 

- 1 > dr. 


( 11 ) 


where Q* is a new calibration constant Including receiver characteristics. 
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The only two unknowns are A (and hence 6R) and O'. Since 1, 

then 






RD 



Hence, If a record could be found In which the radar “saturated", I.e., 
the rain rate was Intense enough to completely attenuate the radar 
returns beyond scmie point, and If A were known, Q' could be calculated. 


VI. CALIBRATION 

A. Matching the Radars 

It was shown In Section V that the calibration problem reduces 
to finding two parameters, A and Q' , for each of the three radars, 
iisiwever, since all the radars are looking at the same rain cell, only 
cn^ pair of values, A and Q' , not necessarily of the same radar, need be 
Kiiown a priori. Using a regression or scatter plot technique, the three 
radars can then be Iteratively bootstrapped Into line and all the other 
cillbratlon constants found. The regression technique used was described 
b> Hodge [7] and Is given below in greater detail. 

B. The Regression Technique 

This technique relies on the principle that when two radars are 
looking at the same event the distributions of the returns of the two 
radars should be the same. As a refinement, returns up to only the point 
of peak reflectivity of the cell are considered. The returns from 
beyond the peak tend to be more seriously corrupted by attenuation 
(Figure Z). 

Using a large number of records, the reflectivity, Zgo, or rain rate, 
R, distributions for the two radars are calculated. These do not, In 
general, coincide (Figure 3). 

Next, the values Z] and Zo in Figure 3 for the two radars are 
computed for several values of % time. Finally, a plot of Z] vs. Zz Is 
made. If the two radars are consistent, the resulting curve should be 
a 45® line. If not, the calibration constants may be adjusted and the 
process repeated till the best possible agreement Is obtained. For 
some results of this procedure see Figures 4, 5 and Reference [7]. 
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Figure 2. Two-frequency data on a typical Figure 3, Equivalent reflectivity factor, 

record. distributions for unmatched rada 
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Figure 4. Comparison between the X and Ku radars obtained using 
the regression-calibration technique. 
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Figure 5. Comparison between the S, X and Ku radars obtained 
using the regression calibration technique. 
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C. Preliminary Calculations 


The equivalent reflectivity factor, Zeq(nw®/m3) , and specific 
attenuation, o (dB/kn), were computed as functions of rain rate, 

R(mm/hr), using Hie theory at all three radar frequencies. The Harshall- 
Palmer drop size distribution and a tei»|)erature of 0“c were assured 
for these calculations. A power law relationship representing these 
values was then developed using logarithmic regression analysis. The 
results are shown in Figures 6 and 7. This process gave the constants 
^z» ®RD* three frequencies. Ttese are given in Table 2. 


Table 2 


Frequency 

(GHz) 

"z 

»z 

®RD ^RO 

3.05 

285.4 

1.453 

9.6x10"^ 0.931 

9.35 

313.1 

1.518 

147.5x10"^ 1.070 

15.56 

326.3 

1.496 

457.2x10“^ 1.063 


D. Obtaining the Radar Constants 

Calculations based on available documentation indicated that the 
gradient, GR, of the analog portion of the Ku band radar receiver up 
to the A/D converter, as it then stood, was approximately 41.8 dBm/volt. 
Thus the A parameter For the Ku receiver, A|(y, was known. 

Next, since the Ku radar would be the first to saturate in heavy 
rain, the value of Q' was calculated for each record of the edited raw 
digital data using Equation (12) and its distribution found. The records 
in which Q’ was highest were then plotted and examined to see if they 
showed evidence of the Ku returns teing completely attenuated by the 
rain beyond soro point, while the S and/or X radar returns still indicated 
a presence of rain past that point. Unfortunately, this was not found, 
indicating that none of the events examined had a sufficiently high 
rain rate or depth to make this happen. Further, the highest value of 
Q* calculated in this manner was still some 15 dB below the theoretical 
value calculated according to Equation (12) using nominal system 
parameters, confirming the above. It should be noted however that this 
does not invalidate the idea that had a sufficiently intense storm been 
recorded, the radar could have been completely calibrated using only the 
A paramter. 

Since this powerful technique for calibrating the Ku band radar 
was not usable with the available data, more heuristic methods were next 
considered. 
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Figure 6. Theoretical attenuation vs. rain rate and the power 
law approximations at 3.05, 9.35 and 15.56 GHz. 
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Figure 7. Theoretical Z 0 q vs. rain rate and power law 
approximations at 3.05, 9.35 and 15.56 GHz. 
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As a first approximation, the known value of A and the best 
available estimate of the Q' value for the Ku radar, based on calcu- 
lations using nominal system parameters, were used and the three radars 
were forced to match consistently using the regression technioue described 
above. It was then found that the Q' value for the S radar (Qg) would 
have to be significantly increased above its theoretical value to match 
the radars properly. Since this implies the unlikely possibility that 
major parameters such as the transmitted power or antenna gain were 
significantly above nwninal values, Qg was next fixed at its nominal 
theoretical value, together with Auy and the regression procedure repeated 
until the radars were again matcheo. This procedure gave values of 
^S» Ax. Qi, QlCy*. other events when the S band radar was not working, 
the Ku radar calibrations were assumed fixed and the X radar matched 
to this (Figures 4 and 5). This completed the calibration procedure. 


VII. ANALYSIS AND RESULTS 

The raw, usable data were converted to the equivalent reflectivity 
factor, Zeq(dBZ), using Equations (10) and (11) and the calibration 
constants described above, and analyzed. Wherever possible, the analysis 
followed that reported earlier by Konrad [8], so that comparisons could 
be made. However, since the rain rates in the usable data were much 
Iwer than those in [8], the rain cells were classified according to the 
peak reflectivity into 8 classes in the 5-45 dBZ range in 5 dBZ intervals. 

The distributions of both the peak reflectivities (Zeg peak) 
also the reflectivities over the whole cell for all the cells were 
calculated for all three radars. These are shown in Figures 8 and 9, 
respectively. Cell width distributions, calculated using fixed 
reflectivity levels, are shown for the three radars in Figure 10. The 
distributions differ because neither of the three radar data bases contains 
the entire data base of the other two radars. The X and Ku radars 
share the most data, however. 

Contours of Zgq at 3 dB intervals below the peak equivalent 
reflectivity were also obtained for each cell. The cell widths measured 
at these contour levels were averaged over each class of cells and are 
shown in Figures 11, 12 and 13. 

The distributions of the reflectivity, Zeq, and the peak reflectivity, 
^eq peak, both show that the cells analyzed did not exceed the 35-40 dBZ 
range, the X and Ku radars show very good agreement. The cell width 
distributions indicate that cells up to 28.5 km wide were observed, though 
infrequently; the three radars agreeing well even over these sizes when 
attenuation would be significant. 

The plots of relative contours vs. cell width clearly indicate that 
as the peak reflectivity increases, i.e., at higher rain rates, the 
dimensions of the cell, at a given reflectivity level below the peak 
decreases. These average approximately 4 km at -6 dB and 11 km at -20 dB 
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Figure 8. Cumulative frequency distributions of peak 
equivalent reflectivity factor. 
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Figure 9. Cumulative frequency distributions of 
equivalent reflectivity factor. 
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Figure 10. Cell width cumulative frequency distributions. 
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Figure 12. Relative reflectivity contour vs. average 
cell width at 9.35 GHz for several cell 
classes. 
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Figure 13. Relative reflectivity contour vs. average 
cell width at 15,56 GHz for several cell 
classes. 
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for the 35-40 dBZ cell class. Tte differences In the dimensions given 
by the three radars are again due, mainly, to the different data bases 
over which these were averaged. 

While It is not possible to make direct comparisons with the results 
reported by Konrad [8] due to the itwch Iwer Intensities observed and the 
limited data, the above results are seen to follow. In general, the trends 
reported by him. Hence, they are consistent with Konrad's results. 

In addition to the above, attempts were made to obtain curves of 
Zpeak vs- cell width and the average slope for each peak reflectivity 
Interval, among other parameters, as reported by Konrad In [8], However, 
due to the limited data base, meaningful results could not be obtained. 


VIII. CONCLUSIONS 

Three- frequency radar test data, taken in 1973, were examined. 
Several methods Including the use of the raw radar data were utilized 
to arrive at reasonable and consistent calibration constants for the three 
radars. Using these constants, cell width, peak reflectivity and 
reflectivity distributions and curves of cell widths at various relative- 
contour levels were obtained. The results are consistent with previously 
reported work. Cell widths vary frwn 4 km at 6 dB below peak reflectivity 
to 11 km at 20 dB below peak reflectivity for the 35-40 dBZ cell category. 
Overall widths of up to 28.5 km were observed. Though the observed rain 
rates did not exceed a corresponding equivalent reflectivity factor of 
40 dBZ, the agreement between the radars was good over even large cells. 

The theory and techniques used follow some new work in this area. 

The results denwnstrate the feasibility of using attenuating frequency 
radars for rain rate measurements. 
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PUBLICATIONS 


In addition to the results In this Final Report, other results of 
the effort under Contract NAS5-23850 were reported in two technical reports. 
The first of these [6] described techniques for radar calibration using 
radiometry and the utility of radar saturation as a nwans of direct, 
absolute radar calibration. The results of this radar calibration analysis 
were also accepted for presentation at the National Radio Science Meeting, 
Boulder, Colorado, January, 1978 [9]. The second technical report [10] 
presented the results of numerical calculations of specific attenuation, 
a, and equivalent reflectivity, ZgQ, for frequencies from 1 to 500 GHz 
and for rain rates from 1.27 to 152.4 mm/hr. The results of a regression 
analysis yielding the constants aj^p, b^^p, a^, and b^ for the pwer law 

relationships, and Zgq=a 2 R^^. were also presented for the same 


frequency range. 
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